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We propose a model for surfaces in mixtures of oil, water, and surfactants with strong 
electric dipoles. The dipole interactions give rise to a non-local interaction with a negative 
stiffness between surface elements. We show that, for large space dimension D, this model 
has a phase transition from rough to smooth surfaces. Contrary to models with a simple 
quadratic curvature energy of a positive sign, which always have a finite persistence length 
typical for rough surfaces, the smooth surfaces in our model exhibit long-range correlations 
of a generalized antiferromagnetic type. These correlations might be related to the recently 
observed "egg carton" superstructure of membranes. 
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Models of random surfaces play an important role in chemical physics, biophysics, and 
particle physics. The statistical mechanics of such models |l|] is usually discussed on the 
basis of a phenomenological Hamiltonian which contains a surface tension and a positive 
extrinsic curvature stiffness 0. Freely suspended vesicles formed from such membranes 
have a vanishing tension, and their shapes are determined by a positive curvature energy 
0. Such a curvature energy is known to have only a finite persistence length ^ 
The correlation functions of the tangent vectors of the surface fall off like e"'^'/^ signaling 
a rough surface on scales larger than ^. For biomembranes, this limitation is of 

no practical importance, since their persistence lengths are much larger than the size of 
the vesicles. This is not true, however, for layers of surfactants between oil and water, 
or for double layers in water alone. Their persistence lengths are of the order of a few 
hundred A. Beyond this distance, membrane fluctuations are governed only by tension. 
This phenomenon is called spontaneous generation of tension 0. Thus, in these models, 
the tension can never disappear completely. Experimentally, however, it is quite easy to 
generate large smooth interfaces. In fact, this is done industrially during the tertiary 
recovery of oil from microemulsions. The fractal surfactant interfaces between oil and 
water are smoothened by the addition of salt. This produces a desired phase separation 
in which the oil floats on top of a single interface. The purpose of this note is to point out 
a mechanism which may explain this. 

Some years ago it was shown that the curvature stiffness of surfactants at the interfaces 
between oil and water [|TD[ receives a significant contribution from the electric fields created 



by layers and bilayers of surface charge [|Tl|. These surfactants possess electric dipoles 
which align along the normals to the interface. The dipoles can be quite large |T^ . 
The range of the electric fields can be varied by the addition of salt. In a field theoretic 
description, the dipoles interact with an electric potential of a finite mass n which is equal to 
the inverse Debye screening length, thus being directly governed by the salt concentration 
n [0: 

where ks is the Boltzmann constant, e the electron charge, and e the dielectric constant. 

In our model we shall ignore the curvature stiffness generated by the mechanical 
properties of the surfactant molecules, or by the van der Waals forces . Thus we consider 
only surfactants with large dipole moments, where the other forces can be ignored, and 
curvature stiffness is completely due to the dipole-dipole interactions. 
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The interaction between the dipole moments in three dimensions is similar to the 
tensor interaction between tangent tensors in the world surfaces of models of quark con- 
finement, in which quarks are held together in the same way as magnetic charges would in 
a superconductor [l^] . Similar models have been conjectured to describe the world surface 
between quarks in the confining phase of QCD ilH]. A recent analysis in Refs. [p!7|,p!8 



indicates that such a model describes smooth surfaces, a fact which seems to be confirmed 
by recent Monte Carlo simulations |T^ . 

By integrating out the electric potential one obtains a non-local interaction between 
normals to the surface. Unfortunetely, such an interaction is hard to treat analytically. 
For this reason we shall study the phase structure of a simplified version of this model 



which retains the essential features of the original one, most notably a non-local interaction 
with a negative stiffness. 

Our starting point is a Hamiltonian describing surfaces with tension, whose normal 
vectors are coupled to a fiuctuating electric field E^{x) = d^4>{x): 

H = r j (f^^ + J^^'^l {d^<l>d^<p + ^'(/>') +dj d^^^ d^4> . (2) 

The surface is parametrized by the three functions x^(^i, ^2) (a* = 1, 2, 3). The parameter 
d denotes the (uniform) dipole density on the surface, and n is the inverse Debye length 
(P. The intrinsic geometry of the surface resides in the induced metric 

gab = daX^dhX^ , 5f = det gab ■ (3) 

The extrinsic geometry is described by the normal vectors 

Na = 77^ ^ab daXf^dbX^ , N^^N^ = 1 . (4) 

Here, greek letters refer to the embedding space, while latin letters denote surface coor- 
dinates. For simplicity we use units in which ^ = 1, c = 1, and measure all momenta, 
energies, or distances in units of the microscopic cutoff A or its inverse, respectively. A 
natural cutoff is provided by the molecular size. We furthermore set A = 1. 

The electric coupling of the surface can be simplified by writing it in terms of the trace 
C = Tr Cab of the second fundamental form, which is equal to twice the mean curvature: 

= r y d^^^ + Jd^x^ {d^<l>d^(j) + +dj d'i^ <j)C . (5) 



Integrating out the electric field, we obtain the non-local energy 

H = rJ d^cVa -yI^^'^I ^ ^"^^^^ " "^^^'^^ Va^Ci^) , (6) 

where 

Y{x) = -—- . (7) 
47r|x| ^ ^ 

is the Yukawa potential. 

We then proceed by rewriting the non-local kernel in the embedding space as a corre- 
sponding non-local kernel on the surface. To this end we introduce a new local coordinate 
system around each point of the surface: ^^,^^,X'^- The coordinates and at = 
describe the original surface. Together with the coordinate x^; they form a locally flat 
coordinate system orthogonal to the surface. The coordinate transformation is described 
by functions x"^) which for x^ = coincide with the original parametrization 

surface. 

The Yukawa potential in can be rewritten as 

JdxJ dx' 5{x) j^^^ X) - x{e, x')) Six') . (8) 

We can now use the transformation rules 

V9 



- = - V 



(9) 



V9 

where Va denotes covariant derivatives along the surface, to rewrite the Yukawa potential 

as 

y{x{o-xm= I dx 6{x) ,\ , s{x)\s\^~o 

where G is formally defined by an expansion in powers of ('D//i)2. In order to compute 
this, we first note that (M^ — V2)~^5(x) is the Yukawa Green's function in one dimension. 



and thus equal to exp (— M|r|) /2M. The first delta function in (|T0[) tells us then that we 
have to take this function at r = 0, implying that 

Inserting this result into (Bf), we obtain 



(12) 



This shows that dipole dominance leads to a non-local interaction with negative stiffness 

The formulation (^2]) is rather awkward for analytic computations due to the non- 
linear character of the trace C. For this reason we shall investigate a model [|2^, which 



is simpler to handle while embodying the most important features of (0), namely non- 
locality and negative stiffness. This model is formulated in terms of the tangent vectors 
VaX^ and posseses the following Hamiltonian: 



/3H^ J (fi^ g'''VaX^W{z,m,V^) Vt 



hX^ , (14) 

Here the index ii runs over ii = 1, . . . , with D the dimension of the embedding space, 
which we keep variable from now on. The interaction kernel is 



W{z,m,V'^) = - ^ where z = ^, m = . (15) 

1 — /m^ 2 a 

It has the expansion 

W {z, m,V^) = z + sV^ + . . . , (16) 

with s = z/m'^ = (3k being the reduced dimensionless stiffness parameter. These expansion 
terms correspond to the Hamiltonian 

= r y d'^v^ +kJ (f^^ g'^^V^Xf.V^V.x^ + ... , (17) 

which matches (up to boundary terms) ( [T^ ) to this order. 
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This model can now be analyzed non-perturbatively by standard large-D techniques 
along the lines of Refs. ||^,^,^ . To this end we introduce a dimensionless Lagrange 
multiplier matrix A"'' to enforce the constraint gab = daX^dtXfj,, 

f3H^(3H + j d'^CVg y'^dax^dbx^-gab) ■ (18) 
We then parametrize the world-sheet in the Gauss map as 

x^{i)={iui2A\0) , i = 3,...,D, (19) 

where -Ri/2 < < Ri/2, -R2/2 < < R2/2 and (p'{^) describe the - 2 transverse 
fluctuations. Then we search for an isotropic saddle point of the form 

9ab = pSab, X"' ^Xg"' , (20) 

for the metric and the Lagrange multiplier of infinite systems {Ri,R2 — > 00). At such a 
saddle point, we obtain a Hamiltonian 

f3H = 2 j (fi [z + A(l -p)] + j (fi daf [X + W (z, m, V^)] da(p' . (21) 

Integrating out the transverse fiuctuations, we obtain in the infinite-area limit the free 
energy 

PF = 2Aext (z + A(l - p)) + Aext^^^P J d^P In {p' [X + W {z, m,p'')] } , (22) 

where A^^t = R1R2 is the extrinsic, physical area. For large D, the fiuctuations of A and 
p are suppressed and these variables take their classical values, determined by the two 
saddle-point equations 

D — 2 

X = / dp plnUp^ \\ + W [^z,m,p'^Y\\ , 

^\ (23) 
p-l D-2 ^ 1 



dp p 



,2^ ' 



p Stt Jq \ + W {z,m,p^) 

where we have introduced the ultraviolet cutoff by restricting the momentum integrations 
to p < 1. Inserting the first saddle-point equation into (|22|) , we find 



(3F = 2{z + X) Aext , 
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(24) 



showing that the tension is renormahzed to 



a = 2{z + X) . 



(25) 



The physics of our surfaces in the large- -D limit is determined thus by the two saddle- 
point equations (^Bp. The first of these equations requires the vanishing of the saddle- 
function 



D-2 

8n 



dp p In {p^ ( A -f [2, m, p^) ] } 



(26) 



and determines the Lagrange multiplier A as the solution of a transcendental equation. 
After this, the second saddle-point equation determines the metric via the equation 
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P 



f'{z,m,X) 

where a prime denotes the derivative with respect to A. 



(27) 



The simple interaction ([T^) in ( |14D has the advantage that the saddle function (|2 
can be computed exactly: 



f{z,m, A) = A 



-1 - m In 1 



D-2 
16% 



D-2 f m^{X + z) + X 



1 



m'^(X + z) ^ fm^iX + z) 
- In 



A 



1 -\- m? 



X 



In A + 



zm 



1 + m? 



(28) 



This function has the following limiting values: 



lim f{z, m, A) = 00 

A — >oo 



lim 

AT-nin 



-zm 



-zm 



D-2 
167r 



(— 1 + In zm^) + O (m^ In m^) , m ^ 1 , (^29) 



1 + m 



2 ■ 



For z sufficiently large we have lim;v^Amin7 < there exists at least one solution to the 
saddle-point equation f{z, e, A) = 0. 

The derivative of the sad die- function /, 



f'{z,e,X) = l 



D-2 1 zm-" 



16tv 



X A2 



m^J A + ^ 



(30) 



determines the metric element p via (p7[). Given that 

hm f' {z, m, A) = 1 , 
Hm f'{z,m,X) = — oo , 

the saddle-function / must have an odd number of extrema. Our numerical analysis shows 
that it has exactly one minimum. When this minimum lies above zero, the saddle-point 
equations have no solutions. When the minimum lies below zero the saddle-point equation 
/(z, m, A) = has two solutions. Only the largest of these two solutions, however, is 
physical since at the smallest one f'{z,m,X) = 1/p < 0. Then we have exactly one 
physical solution of the saddle-point equations. 

In Fig.l we plot the critical line in parameter space below which there exists one 
solution to the saddle-point equations for D = 3. We choose to plot this line as a function 
of the parameters m and ksT/r = l/2z, where T is the temperature. 



0.8- 




4 6 8 10 

m 



Fig. 1: Critical line for roughening transition at large- 



When the critical line is approached from below, the minimum of the saddle-function 
f{z, m. A) goes to zero, so that A ^ A* which is the solution of the equation /' (z, m. A*) = 
0. Using (pTl), we conclude that 

lim p = oo . (32) 
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Thus, when approaching to the critical line, the intrinsic area of the surface diverges 
whereas the renormalized tension remains finite. As far as the phase transition with 
respect to T or 1/r is concerned, this property is also found in the model with positive 
stiffness In the low-temperature phase of that model, however, there are short-range 
correlations between normals to the surface, indicating roughness . As we show below, 
the situation is totally different in our model, as a consequence of the negative stiffness. 
Note that our model reproduces correctly the positive slope of the critical line in (T, 
n)-space, typical of ionic surfactants [jl^ . 

The geometric aspects of the surface are embodied in the correlation function 

9ab{C - a = {da<P\0 d,<P\a) , (33) 
for the normal components of the tangent vectors. This is immediately obtained from (|2l] ) 



as 

9ai>{^ - a = SaB TTT^ / d'p ^ e^Vp^C^'^') , (34) 

(27r)^ J 2[X + W{z,m,p^)\ 

where W (2;, m^p^) is the Fourier transform of the interaction (|15[). This correlation func- 
tion can be rewritten as 

9a,i^ - a =Satj:^ fd'p\^+ ( ' 1,1 X \ ^'^''^'-''^ : (35) 
(27r)2 ] I 2A f^_2s%2^ 



where we have used the fact that the solution for the saddle-point value of A is negative. 
The dominant large-distance behaviour of this correlation function is given by: 



9ah{i - O - Sab 



■ (36) 



zcp \c _ ci\ \\ 2s A 4 



It is easiest to check this result by computing backwards its Fourier transform. To this 
end it is useful to recognize that (|36|) represents the asymptotic behaviour of the von Neu- 
mann function {z/d,s}?)NQ {^^Jzap/2s\\\\^ ~ ^1 j • Computing the two-dimensional Fourier 
transform of this function one finds m ;z/|A|(a - 2s{\\\/z)p^) for all < p < ^za/2s\\\ 
i.e. one reproduces exactly the correlation (|3^) in momentum space. 

The correlation function ( |36| ) is valid up to a large infrared cutoff R such that (1/R) 
provides a regulator for the pole at p = ^J za/2s\\\) in (^5|). So, we have long-range order 
in the sense that there is no finite correlation length and J d'^a 9aa{0 diverges like VR. 
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The oscillatory behaviour is due to the negative stiffness: it indicates that the surface 
fluctuates around the reference plane with an intrinsic wavelength 



^«\/^- (37) 



These long-range correlations are thus of a generalized antiferromagnetic type, with do- 
mains of orientational order of size I. This is a consequence of the frustrated antiferro- 



magnetic interaction with negative stiffness [18 



Contrary to models with a positive stiffness, our model predicts the existence of a 
smooth surface phase with long-range correlations for small T/r and large m. Notice that, 
due to ( ^2|) the correlations (9a0*(O vanish on the critical line.te 

We don't know the exact nature of the roughening transition occuring on the critical 
line depicted in fig. 1. The diverging intrinsic area and the vanishing correlators, however, 
indicate for certain that the rough surface is a fractal object. Note, however, that it is not 
of the type expected for a crumpled surface, where we should have a homogeneous solution 



to the saddle-point equations, while the correlation functions (|3J) would be short-range. 
It is not clear to us whether a generic solution of the saddle-point equations, not restricted 
by the ansatz (pOD, exists beyond the critical line. Such a solution would correspond to a 
fractal phase where the surface develops fingers and bumps. Also, it is not clear whether 
the large- L> transition survives the corrections to D=3. If not, the surfaces would always 
be in the smooth, perturbative phase. In any case, we have shown how dipole dominance 
provides a viable mechanism for large, smooth surfactant interfaces. In addition, we have 
shown how a non-local interaction with negative stiffness leads generically to long-range, 
periodic correlation functions. Thus, our model might be of relevance for the recently 



observed "egg carton" superstructure of membranes, although there are also other 



possible mechanisms to explain these p4 
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